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Along the lines of nonextensive statistical mechanics, based on the entropy Sq = k(l — 
^2iPi)/(l ~ -0 i^i = — k£iPi m Pi)i an d Beck-Cohen superstatistics, we heuristically 
generalize Planck's statistical law for the black-body radiation. The procedure is based 
on the discussion of the differential equation dy/dx = —aiy~{a q —a\)y q (with y(0) = 1), 
whose q = 2 particular case leads to the celebrated law, as originally shown by Planck 
himself in his October 1900 paper. Although the present generalization is mathematically 
simple and elegant, we have unfortunately no physical application of it at the present 
moment. It opens nevertheless the door to a type of approach that might be of some 
interest in more complex, possibly out-of-equilibrium, phenomena. 

We normally obtain the statistical mechanical equilibrium distribution by 
optimizing, under appropriate constraints, an entropic functional, namely the 
Boltzmann-Gibbs (BG) entropy Sbg — —kJ2iPi^ n Pi- The success and elegance of 
this variational method are unquestioned. But at least one more possibility exists, 
namely through differential equations. Such a path is virtually never followed. In- 
deed, such an approach might seem quite bizarre at first sight. But we should by 
no means overlook that it has at least one distinguished predecessor: Planck' s law 
for the black-body radiation. Indeed, Planck published two papers on the subject 
in 1900. The first one in October, the second one in December D The bases of 
both of them were considered at the time as totally heuristic ones, although kind 
of different in nature. The second paper might be considered as a primitive form of 
what has now become the standard approach to statistical mechanics, based on the 
optimization of an entropy functional, the connection with Bose-Einstein statistics, 
and, ultimately, with the Boltzmann-Gibbs thermal theory for a quantum harmonic 
oscillator The first paper "™ however, is totally based on simple arguments re- 
garding an ordinary differential equation. It is along this line that the present paper 
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is constructed. 

If Sbg is extremized under appropriate constraints, we obtain the famous BG 
weight p(E) = p(0) e~P E . This distribution can be seen as the solution of the dif- 
ferential equation dp/ dE = —ftp . Since more than one decade, a lot of effort is 
being dedicated to the study of the so called "nonextensive statistical mechanics" , 
based on the generalized entropy S q = k(l — J2i Pi)/(l ~ 1) (&1 = Sbg) ^ (for a re- 
view, see^f). The extremization of this entropy under appropriate constraints yields 
p(E) = p(0)e-P E , where ef = [1 + (1 - q)x] 1 ^ 1 ~^ (ef = e x ). This distribution, 
which has been shown to emerge in many natural and artificial systems can be 
seen as the solution of the differential equation d[p/p(0)]/dE = — (3[p/p(0)] q . As a 
next step, we may consider even more complex systems, namely those which exhibit, 
for increasing E, a crossover from nonextensive to BG statistics. Such appears to 
be the case of cosmic rays ®. Such situations can be handled with a differential 
equation which unifies the previous two ones, as follows: 

= -AbMo)] - (A, - lh)Wp(P)] q ■ (i) 

Excepting for the fact that here q may be noninteger, this differential equation is a 
particular case of Bernoulli' s differential equation. Its solution is given by 

P(E) = - (2) 

1 + |a ( e (g-l)^B _ 1) 

which precisely exhibits the desired crossover for q > 1 and < /3i « (3 q . Indeed, 
for (q - l)PiE « 1 we have that p/p(0) ~ e q E ', whereas, for (q ~ l)fiiE >> 1, 
we have p oc e~@ lE . 

In the limit f3 q / j3\ — > oo and p(Q)f3i/f3 q — » C, where C is a constant, Eq. (2) 
becomes 

P(E) = - -j-, (3) 

L,(q-l)piE _ I ] «-i 

= ■ ( 4 ) 

If we multiply this statistical weight by the photon density of states g(E) oc E 2 and 
by the energy E, we have the celebrated frequency spectral density 

< v ) a e hu/k B T _ ! ' ( 5 ) 

where we have identified /3i — > 1/fcsT and — > /i^. It is in this precise sense that 
Eq. (3) (hence Eq. (2)) can be seen as a generalization of Planck statistics. 
For q > 1, Eq. (3) can be written as 

^ = f>(n,g) e -^, (6) 

n=0 



which, for q = 2, becomes 



February 2, 2008 12:39 Proceedings Trim Size: 9.75in x 6.5in 



SouzaTsallisErice 



3 




Figure 1. Degeneracy d(n, q) as function of q (n = 0, 1, 2, 3); q = 2 corresponds to Planck law. 



where 



and 



E n = [(q-l)n+l]E <xn- 



d(n, q) = 



l 

9-1 



(^i)r(n- 



1) 



(7) 



(8) 



r(x) being the Gamma function. We may now follow Planck' s path in his December 
1900 paper, where he introduced the discretization of energy that eventually led to 
the formulation of quantum mechanics. Consistently, we may interpret E n as a 
discretized energy and d(n, q) as its degeneracy. We see that, Vg > 1, the spectrum 
is made of equidistant levels, like that of the quantum one- dimensional harmonic 
oscillator. The situation is definitively different in what concerns the degeneracy 
(see Fig. 1). Only for q = 2 we have the remarkable property d(n,2) = 1 (Vn), 
which recovers the harmonic oscillator problem. 

At this point, let us emphasize that any thermostatistical weight (that of ther- 
mal equilibrium for instance) reflects the microscopic dynamics of the system. This 
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fact was addressed by Einstein in 1910 ^ , and was recently revisited by several au- 
thors (see El for instance). It was shown also, on quite general grounds, in^l In 
the same vein, a dynamical theory of weakly coupled harmonic oscillators system 
was recently used for deducing the functional relation between energy variance and 
mean energy that was conjectured by Einstein in connection with Planck' s for- 
mula, thus exhibiting that it is a consequence of pure dynamics It is within this 
dynamical interpretation that Beck and Cohen introduced their superstatisticQ^. 
Indeed, nonequilibrium systems might exhibit spatio-temporal fluctuations of inten- 
sive quantities, e.g., the temperature. They assumed then that the inverse temper- 
ature (3 might itself be a stochastic variable, such that the generalized distribution 
of energy is expressed as 



p(0) 



dpfifi)*-^ , (9) 



U 



where the distribution f(fi) satisfies J °° dj3f{[3) — 1. The effective statistical me- 
chanics of such systems depends on the statistical properties of the fluctuations of 
the temperature and similar intensive quantities. Naturally, if there are no fluc- 
tuations of intensive quantities at all, the system must obey BG distribution (i.e., 
f(0) = S(f3 — 1/fceT)). They also showed that, if /(/3) is the 7-distribution (see 
also'^'), one obtains the ^-exponential weight of nonextensive statistical mechanics. 
Moreover, for small variance of the fluctuations, the nonextensive statistical distri- 
bution is once again reobtained. See'^for an entropic functional which, extremized 
under appropriate constraints, recovers the distribution of superstatistics. 

We straightforwardly obtain, through Laplace transform, that the superstatis- 
tical distribution f((3) corresponding to the p(E)/p(0) given by Eq. (2) is 

f{P) = (f ) 51 S ^ ^ i 1 ~jj 5{f3 ~ MiQ ~ + 1D ' (10) 



Moreover, we define 



(0 2 ) 

where (...) = J °° d(3f (/?)(...). The notation q BC (BC stands for Beck-Cohen) has 
been introduced to avoid confusion with the present q. Only when /(/3) equals the 
7-distribution we have qsc — <?• Using Eq. (10) and integrating we obtain 



<fl)=0 q and tfP)=0$(q-(q-l)^y 
Replacing(12) into (11) we obtain 



(12) 



(13) 
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(a) 




(b) 






(o) 

l 1 1 ■ 






1 1 ■ , ■ 



Figure 2. Functions -^yy (left) and f(/3) (right) for (/3) = 1. (a) Boltzmann-Gibbs distribution 
[£<§> = e -B. m = S(P-l)} ; (b) g = g BC = 1.8 distribution = ^ gfegj ^ 

= 08^rW e_1 ' 25/3 ]; W (9.9BO) = (2,3/2) distribution = /(/?) = 

\[5(P - §) + |5(/3 - 1) + 3<5(/3 - |) + ..-]]; (d) (g,g B c) = (3/2,5/4) distribution = 

1) + ...]J. In the cases (a,c,d), what is 



[4(i-ln2)]- 



/(/?) = \[5(p-\) + 5(p 



I) + |*08 



(2e E / 4 -l) 2 ' 1 V> 4L"V 2' 

represented is not / (/?) strictly speaking, but rather the weights of the Dirac delta's. 



It is worthy remarking that, for all admissible f(0), we can write the asymptotic 
expression p(E)/p(0) = (e"^) ~ e~^ E (l + where a = ^/WT^W = 

Finally, we may rewrite distribution (2) as follows: 



P(E) = L g-1 
hence, through Laplace transform, 



,(q-q BC )(l3)E^ 



(14) 



/(/?) = 



1-9 



^d(n,gr) 



n=0 



(9-1) 



[(g-l)n+l] 



(15) 

Observe that, for all q, if q BC 1 we obtain the BG distribution. In addition, 
we see that j3 generically assumes discrete values in /(/?) . If we focus on the limit 
of continuous values for (3, we must have (using Eq. (10)) A/3 = (3(n + 1) — f3(n) = 
(3i(q — 1) — > 0, and this is obtained (see Eq. (13)) when (/3) — ► (i.e., high 



2, 2008 12:39 Proceedings Trim Size: 9.75in x 6.5in 



SouzaTsallisErice 



6 

temperature) or q BC — » q (i.e., g-statistics) . In Fig. 2 we present typical examples 
of pairs (p(E)/p(0), /(/?)). 

Summarizing, we obtained the distribution corresponding to the differential 
equation (1), expected to characterize a class of physical stationary states where a 
crossover occurs between nonextensive and BG statistics. This led us to a possi- 
ble generalization of Planck law. We obtained also the Beck-Cohen superstatistical 
distribution /(/?) associated with such type of crossovers between statistics. Along 
similar lines, it is possible to study crossovers between q and q' statistics, with 
eventual applications in turbulence and other complex phenomena. 
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